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ABSTRACT 


Two  experimental  techniques  to  determine  the  frequency- 
dependent  complex  mass  of  a  viscous  fluid  contained  in  a  rigid  porous 
solid  are  investigated.  In  one  technique  the  moving  mass  of  the  fluid 
is  sensed  by  its  effect  on  the  measured  input  electrical  impedance  of  a 
moving  coil  transducer.  In  the  second  technique  the  moving  mass  is 
extracted  from  the  measured  acoustic  pressure  required  for  the  fluid 
to  oscillate  with  a  known  amplitude  through  the  solid  frame.  A 
description  of  the  apparatus  and  preliminary  results  using  the 
pressure  method  as  well  as  a  comparison  of  the  two  techniques  are 
presented. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Acoustic  propagation  in  the  ocean  is  frequently  controlled  by  the 
sea  bottom  (Urick,  1982,  p.  11-1).  Reflective  bottom  loss  due  to 
impedance  mismatch,  scattering  due  to  bottom  roughness,  and 
attenuation  due  to  dilitational  and  shear  wave  excitation  in  bottom 
sediment  all  dominate  shallow  water  transmission  at  low  frequencies, 
and  influence  deep  water  propagation  to  short  and  intermediate 
ranges  at  mid-frequencies  (Tolstoy  and  Clay,  1972,  p.  8).  The 
importance  of  these  physical  phenomena  to  underwater  acoustics  is 

# 

well  recognized.  Extensive  literature  is  available  on  the  broad  subject 
^  of  marine  sediment  acoustics  because  of  the  vast  contributions  by 

many  renowned  individuals  over  the  past  several  decades,  e.g., 
Hampton  (1974),  Hamilton  (1974),  Kuperman  and  Jensen  (1980), 
Akal  and  Berkson  (1986),  and  Stoll  (1986). 

Even  so,  one  example  of  bottom  interaction  physics  that  remains 
difficult  to  specify  theoretically  and  validate  experimentally  is  the 
frequency  dependence  of  the  attenuation.  Because  of  the  difficulty  of 
measuring  attenuation  at  low  frequencies,  values  are  obtained  by 
extrapolation  from  accurate  high  frequency  measurements;  too  often 
decades  removed  from  the  frequency  span  of  study  (Hovem,  1980,  p. 
2).  This  problem,  coupled  with  the  inherent  experimental 
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uncertainties  and  the  fact  that  no  single  data  set  spans  the  entire 
frequency  range  of  interest  in  underwater  acoustics,  i.e.,  Hz  to  MHz 
(Kibblewhite,  1989,  p.  718),  underscores  the  requirement  for  an 
accurate  theoretical  model  for  the  attenuation  of  sound  in  sediments. 
Such  a  theoretical  model  can  be  prescribed  by  applying  the  Biot  theory 
of  acoustic  wave  propagation.  Biot  theory  accounts  for  fluid-solid 
viscous  interaction,  which  is  an  important  ingredient  in  the 
phenomena  of  bottom  attenuation.  However,  the  details  of  this 
interaction  are  not  well  characterized  for  porous  solid  geometries. 

The  Biot  model  of  fluid-saturated  sediment  contains  thirteen 
parameters,  defined  in  Biot  (1956a,  1956b),  Stoll  (1974,1985,1986), 
and  Plona  and  Johnson  (1984),  which  depend  on  the  microscopic 
pore  geometry  of  the  porous  solid,  the  bulk  properties  of  the  solid, 
and  the  bulk  properties  of  the  fluid.  These  parameters  are  very 
difficult  to  extract  from  conventional  measurements  of  compressional 
sound  speed  and  attenuation.  This  fact,  along  with  the  paucity  of 
independent  observations  of  the  effects  of  fluid-solid  interaction,  has 
been  one  of  the  major  obstacles  to  general  acceptance  of  the  Biot-Stoll 
view  by  the  underwater  acoustics  community  (Kibblewhite,  1989,  p. 
718).  However,  recent  developments  include  success  in  defining 
some  of  the  key  Biot  parameters  and  coefficients.  As  reported  by 
Kibblewhite  (1989);  Bedford  et  al  (1984),  Johnson  et  al  (1987),  and 
Baker  (1986)  have  contributed  significantly  to  specifying  parameters 
required  as  input  to  the  Biot  formulation  of  elastic  wave  propagation  in 
a  fluid  saturated  porous  solid. 
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Least  well  known  are  the  three  parameters  which  depend  on  the 
microscopic  pore  geometry  and  which,  along  with  the  bulk  fluid 
density,  completely  determine  the  so-called  complex  mass  density,  P  , 
of  the  fluid  in  the  pores  (Ogushwitz,  1985,  p.430,  and  Stoll,  1985, 
p.  1798)(as  revealed  in  Chapter  II,  the  complex  mass  density  for  a 
rigid  and  stationary  porous  solid  is  the  negative  of  the  ratio  of  the 
gradient  in  fluid  pressure  to  its  acceleration).  For  this  reason,  a 
measurement  of  P  is  a  superior  probe  of  these  parameters  (Baker, 
1986,  p.  2). 


B.  OBJECTIVES 

The  goal  of  this  research  program  is  to  determine  the  frequency- 
dependent  complex  mass  of  a  viscous  fluid  in  a  rigid  and  stationary 
porous  solid  for  various  typical  pore  geometries.  Two  experimental 
techniques  are  being  investigated,  with  the  fluid-saturated  solid 
sample  contained  in  a  small  cylindrical  cavity,  depicted  in  Fig.  1 
(Grant,  1988,  p.  4).  The  first  technique,  the  so-called  impedance 
method,  is  based  on  electroacoustic  network  theory.  The  moving 
mass  of  the  fluid  is  sensed  by  the  measured  input  electrical  impedance 
of  moving  coil  drivers.  Preliminary  results  for  the  impedance  method 
have  been  reported  by  LCDR  Grant  (Grant,  1988,  p.  57).  The  second 
technique,  the  so-called  direct  method,  is  based  on  Newton's  Second 
Law  for  a  viscous  fluid  contained  in  a  rigid  porous  solid.  Stiff 
piezoceramic  drivers  are  used  to  apply  a  known  displacement  for  a 
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given  voltage  input.  The  moving  mass  is  extracted  from  the  observed 
pressure  difference  across  the  sample  required  for  the  pore  fluid  to 
oscillate  with  a  known  amplitude. 


■ 

III 

Figure  1.  Simple  Model. 

The  specific  objectives  of  the  research  described  in  this  thesis  are: 
(1)  to  develop  and  calibrate  an  apparatus  to  extract  the  complex  fluid 
mass  contained  in  a  rigid  porous  solid  from  acoustic  pressure  and 
particle  velocity  measurements  (direct  method),  (2)  to  compare  these 
calibration  results  with  those  of  the  impedance  method  (Grant,  1988, 
pp.  44-56),  and  (3)  to  identify  the  better  method  of  the  two  for 
further  investigation  with  various  porous  solid  sample  geometries. 
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II.  THEORY 


A.  BIOT-STOLL  THEORY  OF  SOUND  PROPAGATION  IN  A  FLUID- 

FILLED  POROUS  SOLID 

As  summarized  by  Stoll  (1974),  propagation  of  sound  in  the  ocean 
bottom  is  mainly  controlled  by  two  fundamentally  different 
mechanisms:  the  first  is  the  friction  at  grain-grain  contacts  due  to  the 
inelasticity  of  the  solid  frame,  and  the  second  is  the  viscous  drag 
between  the  pore  fluid  and  solid  frame.  Sediments  of  low 
permeability  also  exhibit  another  kind  of  viscous  damping,  attributed 
to  local  fluid  motion  around  intergranular  contacts  (Stoll,  1985,  p. 
1790). 

In  principle,  the  motion  of  solid  and  fluid  components  can  be 
described  by  differential  equations  for  a  microscopic  volume  element 
of  porous  material,  subject  to  the  appropriate  boundary  conditions  at 
the  solid  frame-pore  fluid  interface.  Solutions  to  these  equations  can 
be  found  only  if  the  microscopic  pore  geometry  is  completely 
specified,  a  practical  impossibility.  However,  a  phenomenological, 
macroscopic  (averaged  over  the  microscopic  structure)  description  of 
motion  is  possible,  as  long  as  the  acoustic  wavelength  is  much  greater 
than  a  characteristic  pore  radius.  Biot  derived  just  such  a  description 
of  motion  in  his  series  of  classic  papers  (Biot,  1956a,  1956b,  1962). 
Many  equivalent  formulations  of  Biot  equations  and  parameters  exist  to 


5 


suit  particular  applications.  My  choice  is  the  set  of  equations  of 
motion  presented  by  Baker  (1986,  p.  29): 


P U  =  -P  Au  -u)  +  bv(v-u)  +  cv(v  d  ) 

dt2  dt2 

-  NVxVxu  -  bF^(u  -  U) 

(2.1) 


pj^Us-pAu  -u)  +  CV(V  u)+DV(V  U)  -bF|(U-u) 
dt2  dt2  <*  (2.2) 


where. 


U  =  average  fluid  displacement, 
u  =  average  frame  displacement, 


P  i  =  mass  density  of  bulk  solid  multiplied  by  the  volume 
fraction  of  solid,  i.e.,  Pi  =(1  -P)ps , 

P2  =  mass  density  of  bulk  fluid  multiplied  by  the  volume 
fraction  of  fluid,  i.e.,  P2  =  PPf , 

Pa  =  effective  mass  density  normalized  to  the  apparent 
fluid  density,  i.e..  Pa  =  (a  -  l)Ppf , 

P  =  porosity,  defined  as  the  volume  fraction  of  fluid, 

a  =  tortuosity  parameter,  a  function  of  the  microscopic 
pore  geometry, 

b  =  friction  coefficient  =  PR^c  ;  Rdc  is  the  DC  flow 

resistance,  defined  as  the  average  gradient  in  fluid 
pressure  divided  by  the  average  fluid  velocity  in  the 

direction  of  -^P, 

N  =  shear  modulus  of  the  frame, 
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F  =  complex,  frequency-dependent,  universal  correction 
function. 


The  coefficients  B,  C,  and  D  are  generalized  elastic  constants,  and  are 
functions  of  the  bulk  fluid  and  solid  elastic  contants,  and  the  porosity, 
as  follows: 


B  = 


l-P-!4l-PJK,  +  P^Kb 

K*J  *Sf _ +  4-n 


1P 


+p! 


(2.3) 


PKS 

1  .p.^b 

+  P§, 

Ks. 

Kf 

P2KS 

1  -p.^b 

+  P^s 

KSJ 

Kf 

where, 


(2.4) 


(2.5) 


Ks  =  solid  grain  bulk  modulus, 

Kb  =  solid  frame  (absent  the  pore  fluid)  bulk  modulus. 


Kf  =  pore  fluid  bulk  modulus. 

For  the  case  of  a  stationary  frame,  to  which  this  study  is  restricted, 
u  =  0.  Noting  that  Pa  +  P2  =  «Ppf ,  Eq.  (2.2)  yields,  for  u  =  0: 


aPpf^-U  =  DV(V  •  U )  -  bF— 

at2  at  (2  6) 

For  a  perfectly  rigid  frame  Kf «  Ks,Kb  t  and  Eq.  (2.5)  reduces  to 
D  =  PKf  .  Recalling  also  that  b  =  PRoc  ,  Eq.  (2.6)  becomes: 
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(2.7) 


aPpAj  =  PKfV(  V  •  U )  -  PRdcF^- 
dt2  <* 

Assuming  harmonic  time  dependence  (ei®*)  ,  Eq.  (2.7)  becomes: 

-CD^apf  +  ^^f|U  =  KfV(V  •  U  ) 

From  the  linearized  continuity  equation. 


.dU 


V-U 

pf 

Applying  this  result  and  the  definition 


Kf  =  PfJ* 

8pf 

for  the  bulk  modulus  to  Eq.  (2.8)  gives: 


(2.8) 


(2.9) 


(2.10) 


co^apf  +  =  *^P 


(2.11) 


Defining  the  complex  mass  density  by 


p  =  apf  +  R:^F 


(2.12) 


Eq.  (2.9)  can  be  cast  in  terms  of  fluid  particle  velocity  u  =  jcoU  as 

jcopu  =  -Vp  (2.13) 

which  is  in  the  familiar  form  of  Euler's  equation  (KFCS,  1982,  p.  104), 
(hereafter  the  subscript  f  on  P  will  be  omitted:  it  is  to  be  understood). 
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The  real  part  of  Eq.  (2.12)  is  the  apparent  fluid  mass  density  and 
the  imaginary  part  is  the  negative  of  the  apparent  flow  resistance 
divided  by  the  angular  frequency: 


Mp)  =  Prflfto)  *  R«(“pf T  ^f) 

(2.14) 

ta(?)=1R^M  =  InlJap(+toF| 

(2.15) 

Both  Eq.  (2.14)  and  Eq.  (2.15)  depend  on  frequency.  The  normalized 
complex  mass  density,  p/Pf  ,  is  called  the  complex,  or  dynamic, 
tortuosity,  a  .  and  depends  only  on  the  microscopic  pore  geometry 
and  frequency.  In  the  context  of  the  Biot  theory,  the  microscopic 
pore  geometry  is  characterized  by  the  three  parameters:  (intrinsic) 
tortuosity  a,  DC  flow  resistance  Rpc.  and  the  structure  factor  5,  which 
is  a  scale  factor  for  the  argument  of  the  complex  universal  correction 
function  F.  In  Biot's  original  notation,  these  correspond  to  the  pore 
size  parameter,  permeability,  and  structure  factor,  respectively 
(Ogushwitz,  1985,  p.  430). 

Viscosity  of  the  pore  fluid  is  very  important  to  solid/fluid  coupling. 
A  pertinent  quantity,  therefore,  is  the  viscous  penetration  depth  (5v), 
or  skin  depth,  defined  as 


(2.16) 


9 


where  rj  is  fluid  (shear)  viscosity,  (e.g.,  5V  =  56p.,  for  water  at  20°C  and 
<o  =  100  Hz).  The  relative  magnitude  of  the  imaginary  and  real  parts  of 
the  complex  mass  density,  and  therefore  the  relative  importance  of 
fluid  viscosity  versus  inertia  in  determining  the  fluid  motion,  depends 
upon  the  ratio  of  the  viscous  penetration  depth  to  a  characteristic 
pore  size.  If  5V  is  small  compared  to  a  characteristic  pore  size,  then 
the  motion  of  the  fluid  is  very  similar  to  that  of  an  inviscid  fluid,  and 
independent  (of  the  solid)  sound  propagation  within  the  fluid  is 
possible.  This  case  is  termed  the  high  frequency  limit.  By  contrast,  if 
the  viscous  penetration  depth  is  large  compared  to  the  characteristic 
pore  size,  the  fluid  is  essentially  locked  to  the  rigid  solid  frame.  In 
this  case,  sound  cannot  propagate  independently  within  the  fluid,  and 
it  is  referred  to  as  the  low  frequency  limit.  (Bourbie  et  al,  1987,  p. 
86) 

Two  methods  for  experimentally  determining  the  complex  mass 
density  of  a  fluid- saturated  porous  solid  are  being  investigated,  termed 
the  direct  method  and  the  impedance  method.  A  theoretical 
description  of  each  follows. 


B.  DIRECT  METHOD 

In  the  so-called  direct,  or  Newton's  Second  Law,  method,  the 
complex  mass  density  p  is  directly  obtained  from  a  measurement  of 
the  acoustic  pressure  gradient  required  to  produce  a  certain  (known) 
acceleration,  using  Eq.  (2.13): 
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~  -Vp 
p  =  — ^ 

J®u  .  (2.17) 

The  description  of  how  pressure  and  acceleration  are  measured  is 
given  in  the  next  chapter. 


C.  IMPEDANCE  METHOD 

In  the  impedance  method,  the  complex  mass  density  is  indirectly 
obtained  from  a  measurement  of  the  input  electrical  impedance  (Zjn)  of 
a  moving  coil  transducer  as  follows,  from  Wilson  (1985,  p.  28): 

Zin  =  Ze  + - - - 

Rm  +  Rr  +  jcajM+Mf)  +  78- 

J®  ,  (2.18) 

where: 

B1  =  transduction  coefficient  of  the  driver, 

Rm  =  mechanical  resistance  of  the  driver, 

Rr  =  effective  mechanical  resistance  of  the  fluid, 

M  =  moving  mass  of  the  driver, 

H  =  effective  mass  of  the  fluid, 
s  =  stiffness  of  the  driver, 

Ze  =  blocked  electrical  impedance  of  the  driver, 

The  effective  mass,  Mr  ,  and  the  effective  mechanical  resistance,  Rr  ,of 
the  fluid,  obtained  from  the  imaginary  and  real  parts  of  Eq.  (2.18), 
respectively,  are  given  by: 

1 1 


(2.19) 


M,  =  £lm 


(BlP 

Zjn-Ze. 


+  -5--M 

CD2 


and 


Rr=  Re 


(Blp 


Zis-Ze 


-R, 


m 


(2.20) 


Finally,  H  and  Rr  are  related  to  and  Rfiow  ,  respectively,  by 


„  _  R  _  Rr 

Peff-PSL  and  Rflow  PSL  ,  (2.21) 

where; 

P  =  the  porosity, 

S  =  the  cross-sectional  area  of  the  chamber, 

L  as  the  length  of  the  fluid  plug, 

PSL  =  the  volume  occupied  by  the  fluid. 


A  description  of  the  apparatus  developed  to  extract  the  complex  mass 
from  impedance  measurements  using  a  moving  coil  loudspeaker 
follows  in  the  next  chapter. 
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m.  EXPERIMENTAL  DESIGN.  PROCEDURE.  AND  PREL1 

RESULTS 


•  i»C 


[ABX 


A.  BACKGROUND  AND  PREVIOUS  WORK 

In  Chapter  II,  theoretical  relationships  were  derived  between 
physically  measurable  quantities  by  which  the  complex  moving  mass  of 
a  fluid  contained  in  a  rigid  porous  solid  can  be  extracted.  This 
chapter  describes  the  design  and  preliminary  results  of  the  two 
experimental  techniques  being  investigated  for  this  purpose. 

For  his  Master's  research,  LCDR  Grant  developed  the  apparatus  for 
the  impedance  method,  (Grant,  1988),  and  his  work  is  summarized 
here.  The  development  of  the  impedance  method  was  accomplished 
in  two  phases:  (1)  selection  and  characterization  of  two  matched 
moving-coil  transducers,  and  (2)  design  and  construction  of  a  test 
chamber.  A  pair  of  Philips  model  AD062T8  one-inch  diameter 
tweeters,  depicted  in  cross-section  in  Fig.  2  (Grant,  1988,  p.  18), 
were  used  as  drivers.  To  ensure  oscillatoiy  plug  flow  of  the  fluid 
through  the  sample  of  porous  material,  the  pair  of  transducers  were 
driven  in  parallel  and  out-of-phase  with  respect  to  each  other  in  a 
system  closed  to  ambient  pressure.  To  satisfy  the  rigid  frame 
condition,  a  massive  six  inch  outer  diameter  brass  test  chamber,  with 
a  centered  1.5  inch  diameter  cylindrical  cavity,  was  designed  and  built 
as  shown  in  Fig.  3  (Grant,  1988,  p.  40).  Porous  samples  are  to  be 
housed  within  a  cylindrical  metal  sleeve,  and  the  sleeve  inserted  into 


13 


Surround 


Dome 


Voice  Coil 

Figure  2.  Philips  Driver. 


the  cavity.  Overall  test  chamber  length  was  limited  to  three  inches  to 
avoid  acoustic  resonances  in  the  experimental  frequency  range  of 
interest,  i.e.,  below  one  kilohertz. 


Vent 


Connection 


Fill 

Connection 


Figure  3,  Cutaway  View  of  Impedance  Method  Apparatus. 
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In  a  series  of  experiments,  with  the  objective  to  calibrate  the 
apparatus  with  the  sample  cavity  filled  soley  with  water,  LCDR  Grant 
successfully  demonstrated  that  the  moving  mass  load  presented  to  the 
face  of  one  moving  coil  transducer  could  be  extracted  from  Its 
measured  input  electrical  impedance.  This  required  a  meticulous  and 
thorough  investigation  of  the  five  properties  of  the  moving  coil  driver 
which  appear  in  Eq.  (2.18)  (the  details  of  which  are  presented  in 
Grant,  1988).  Reported  experimental  uncertainty  in  all  five  quantities 
was  consistently  less  than  4.5  per  cent  (Grant,  1988,  pp.  21-33). 
High-precision  input  electrical  impedance  measurements  were  made 
using  a  Hewlett-Packard  Low  Frequency  Impedance  Analyzer  Model 
4192  (Hewlett-Packard  Corporation,  San  Jose,  Ca.).  LCDR  Grant 
concluded  that  it  appeared  possible  to  extract  the  effective  mass  with 
an  accuracy  of  about  five  percent  over  a  significant  portion  of  the 
frequency  band  100  to  1000  Hz;  e.g.,  Fig.  4  (Grant,  1988,  unpublished 
results). 
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Radiation  Mass  for  20  ml  of  Added  Water 


Figure  4.  Impedance  Method  Mass  Extraction  for  20  ml  of  Added  Water. 


B.  FURTHER  IMPEDANCE  METHOD  EXPERIMENTS 

LCDR  Grant’s  results  using  the  impedance  method  revealed  a 
systematic  error  with  increasing  frequency  in  the  band  of  interest,  as 
Fig.  4  indicates.  The  magnitude  of  this  error  was  found  to  be  greater 
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for  larger  added  mass.  A  considerable  effort  was  therefore  made  In 
the  present  investigation  to  improve  the  Impedance  technique.  A 
discussion  of  the  source  of  the  error  and  a  description  of  the  efforts 
made  to  correct  it  follow. 

Fig.  5  (Grant,  1988,  p.  15)  and  Fig.  6  (corrected  figure  from  Grant, 
1988,  p.  14)  show  the  block  diagram  and  lumped  element 
representation  of  the  equivalent  electrical  circuit  for  the  impedance 
method.  The  subscripts  M  and  R  indicate  the  electrical  equivalents  of 
the  transducer  open  circuit  mechanical  impedance  and  the  radiation 
impedance,  respectively,  in  the  mobility  analog  (KFCS,  1982,  p.  360). 
In  this  analogy,  a  greater  mechanical  or  radiation  impedance  results  in 
a  smaller  equivalent  electrical  impedance.  For  frequencies  below  one 
kilohertz,  the  impedance  of  the  inductance  (representing  mechanical 
compliance)  is  greater  than  any  other  element,  and  this  component  is 
largely  ignorable.  Hence,  for  increasing  frequency,  the  radiation  and 
mechanical  impedances  appear  more  and  more  as  short  circuits. 
Under  these  conditions,  the  degree  to  which  the  blocked  electrical 
impedance  is  known  determines  the  precision  with  which  the 
radiation  impedance  can  be  extracted  from  the  measured  input 
electrical  impedance. 

LCDR  Grant  realized  the  blocked  electrical  impedance  cannot  be 
represented  as  a  simple  resistor  and  inductor  in  series  (as  in  Fig.  6). 
In  fact,  a  frequency  dependent  blocked  electrical  impedance  was  used 
to  obtain  the  results  shown  in  Fig.  4.  Because  it  is  virtually  impossible 
to  completely  block  a  transducer  without  damaging  it,  however,  the 
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value  of  the  blocked  electrical  impedance  which  Grant  used  in  his 
analysis  was  obtained  from  electrical  impedance  measurements  on  a 


Figure  5.  Equivalent  Circuit. 


Figure  6.  Lumper  Element  Equivalent  Circuit  -  Mobility  Analogy. 

similar  moving  coil  transducer  which  had  been  completely  potted  with 
epoxy  (Grant,  1988,  p.  22).  For  the  extraction  of  moving  masses  of  50 
g  or  so,  the  blocked  electrical  impedance  must  be  known  with  much 
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more  certainty,  i.e.,  for  the  particular  driver  being  used. 
Furthermore,  because  the  mechanical  resonance  of  the  driver  with  no 
added  mass  occurs  at  about  1000  Hz,  some  of  the  error  in  the 
extracted  value  for  the  larger  masses  could  be  due  to  insufficient 
knowledge  of  the  mechanical  impedance.  In  summary,  it  is  believed 
that  extension  of  the  impedance  method  to  successfully  extract 
masses  of  50  g  or  more  hinges  on  knowing  the  transduction 
coefficient,  mechanical  resistance,  moving  mass  of  the  driver,  stiffness 
of  the  driver,  and  the  blocked  electrical  impedance,  to  a  higher 
degree  of  accuracy. 

The  first  effort  to  improve  the  impedance  method  involved  finding 
the  best  values  of  these  parameters,  that  is,  the  values  which  gave  the 
best  results  for  the  extracted  mass  for  a  collection  of  known  added 
masses.  A  computer  program  was  written  to  control  an  HP  4192  to 
measure  the  input  electrical  impedance.  Impedance  measurements 
were  made  for  added  masses  of  0  to  100  g  in  10  g  increments  from  10 
to  2000  Hz  in  10  Hz  steps.  A  grid -search  least  squares  program 
adapted  from  Bevington  (1969,  pp.  208-213)  was  used  to  find  values 
of  the  "known"  transducer  parameters  (Bl,  Zm.  Ze)  which  minimized 
the  magnitude  of  the  fractional  error  in  the  computed  and  actual 
radiation  impedance  (jaM,), 


(Bl*2 


(Zin  -  Ze) 


"  Zm 


j0)Mr 


(3.1) 
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and  to  store  the  results  in  a  file  for  use  by  graphics  routines.  Seed 
values  for  the  parameter  values  were  taken  from  Grant  (1988,  pp. 
13,39).  The  results  showed  no  significant  improvement  over  Grant's 
results. 

As  an  alternative  attempt  to  improve  the  impedance  method,  the 
problem  of  finding  the  optimum  values  of  the  transducer  parameters 
in  Eq.  (2.18)  was  re-cast  so  as  to  find  the  complex,  frequency 
dependent  transduction  coefficient,  Bl,  that  gives  the  least  variance  of 
the  mechanical  impedance  of  the  transducer  extracted  from 
impedance  measurements.  From  Eq.  (2.18),  the  mechanical 
impedance  may  be  obtained  from  measurements  of  input  electrical 
impedance  with  a  mass  added,  Zjn  ,  and  with  zero  mass  added, 
Zin.no  mass  ,  according  to: 


Zn,  = 


Z  IN,  NO  MASS  -  Zin 


(Bl) 


(Zm  +  jCoM,) 


-1 


which  has  the  root 


(3.2) 


7  _  jcoMr  , 
Zm - 2 


(jcoM,) : 


i  L 
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(3.3) 


4  (Z  DM,  NO  MASS  -  Zin)J 
This  second  attempt  to  obtain  better  values  for  the  transducer 
parameters  produced  results  within  the  experimental  error  reported 
by  Grant  (1988,  p.  30)  in  the  band  500  to  1500  Hz.  For  frequencies 
greater  than  1500  Hz,  however,  the  optimum  value  of  Bl  was 
unreasonably  low,  for  reasons  unknown.  Following  the  disappointing 
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results  of  attempts  to  improve  the  impedance  method  using  the 
available  apparatus,  it  was  concluded  that  better  results  were  not  likely 
to  be  obtained  without  changing  to  another  driver.  Further 
investigation  focused  on  the  development  of  the  direct  method, 
described  next. 


C.  DIRECT  METHOD  DESIGN 
1.  Closed  System  Design 

Figure  7  shows  a  cutaway  view  of  the  new  apparatus  designed 
and  built  to  extract  the  complex  fluid  mass  from  acoustic  pressure  and 
velocity  measurements  (direct  method).  The  sample  chamber  is 
unchanged  from  the  impedance  method  apparatus.  The  direct 
method  requires  stiff  transducers,  driven  in  such  a  way  as  to  cause  the 
fluid  in  the  sample  to  oscillate  as  a  plug.  For  this  reason,  a  "stack"  of 
piezoceramic  rings  was  chosen  for  the  motor  element  of  each  driver. 
Axially  polarized  Navy  Type  I  lead  zirconate-titanate  rings,  .75  inches 
in  diameter  and  .25  inches  thick,  with  a  .25  inch  diameter  center 
hole,  were  procured  (Channel  Industries,  Inc.,  Santa  Barbara,  Ca.). 
The  piezoelectric  strain  constant  in  the  direction  of  polarization,  d33  , 
was  measured  for  each  element  using  a  Berlincourt  d33  meter  (Channel 
Products,  Inc.,  Chesterland,  Oh.).  From  four  selected  d33-matched 
pairs,  two  nearly  identical  piezoceramic  transducer  stacks  were 
assembled  by  cementing  together  four  rings  for  each  driver. 
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A  piston  and  surrounding  diaphragm  assembly  seals  the  cavity  at 
either  end  of  the  sample.  An  Impedance  head,  B&K  Model  8001 
(Bruel  &  Kjser,  Foster  City,  Ca.),  was  positioned  between  the 
piezoceramic  stack  and  this  assembly.  An  impedance  head  is  a 

v*nt  Piezoceramic  "stack" 


Fill 

Connection 


Figure  7.  Cutaway  View  of  Direct  Method  Apparatus. 

combination  of  a  force  gauge  and  an  accelerometer  in  one  housing. 
This  precision  device  allowed  simultaneous  measurement  of 
acceleration  and  force  very  close  to  the  driving  point  of  the  assembly. 
The  impedance  head  also  secured  the  threaded  end  of  the 
compression  bolt,  providing  a  static  pre-stress  on  the  projector  stack. 
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The  opposite  end  of  each  stack  was  cemented  to  a  solid  brass  slab  five 
inches  in  diameter  and  one  inch  thick. 

2.  Closed  System  Acceleration  Calibration 

The  following  procedure  was  developed  to  match  the 
displacements  of  the  piezoelectric  drivers,  to  ensure  plug  flow  of  the 
fluid  in  the  chamber.  The  experimental  set-up  used  is  depicted  in 
Fig.  8.  An  approximately  constant  voltage  was  simultaneously  applied 
to  each  stack  and  the  resulting  accelerations  were  compared.  This 
was  done  from  10  to  2000  Hz  in  ten-Hz  frequency  steps.  At  each 
frequency,  the  acceleration  of  one  stack,  as  measured  from  the 
accelerometer  in  its  impedance  head,  was  adjusted  to  match  the  other 
stack's  acceleration  in  magnitude  (within  +/-  2%)  by  manually  varying 
the  power  amplifier.  Phase  adjustment  was  not  required.  A 
computer  program  (code  listing  appears  in  Appendix  B)  was  written 
specifically  to  accomplish  the  tasks  of:  (1)  measuring  and  comparing 
the  acceleration  magnitude  of  the  two  stacks  at  each  frequency,  (2) 
prompting  the  operator,  when  necessary,  to  adjust  one  stack  to 
achieve  a  match  within  the  desired  tolerance,  (3)  recording  the 
acceleration  and  force  for  each  stack,  and  (4)  stepping  up  the 
frequency  to  repeat  the  process.  An  HP  9000  Series  300  Computer 
commanded  the  HP  3421 A  Data  Aquisition/Control  Unit  (DAC)  and 
triggered  the  HP  4192A  Low  Frequency  Impedance  Analyzer.  With  a 
complete  data  base  of  force  and  acceleration  magnitude  and  phase  as  a 
function  of  frequency,  the  displacement  of  the  transducers  was 
considered  matched  and  known  to  within  one  percent. 
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Figure  8.  Transducer  Calibration  Experiment  Set-up. 


D.  DIRECT  METHOD  EXPERIMENTS 
1.  Introduction 

Several  preliminary  experiments  were  conducted  using  a  fluid 
filled  sample  chamber  with  the  objective  to  calibrate  the  direct 


method  apparatus.  These  experiments  were  performed  with  a  single 
driver  coupled  to  the  test  chamber  In  an  "open  system"  configuration, 
as  illustrated  in  Fig.  9.  In  Appendix  A,  it  is  shown  that  the  radiation 
load  presented  to  the  open  end  is  neglible  in  comparison  to  the 
moving  mass  load.  Thus,  no  correction  needs  to  be  made  to  the 
theory  in  Chapter  II  for  the  open  system.  With  the  characteristics  of 
the  piezoelectric  drivers  established,  a  suitable  means  to  measure  the 
acoustic  pressure  resulting  from  the  oscillating  drive  was  developed. 
Finally,  a  set  of  calibration  experiments  was  performed,  the  results  of 
which  were  used  to  extract  the  moving  mass  of  several  "unknown" 
fluid  loadings. 

2.  Measurement  of  Acoustic  Pressure 

Referring  to  Fig.  7,  the  original  intention  was  to  measure  the 
acoustic  force  on  the  face  of  each  driver  using  the  impedance  head. 
However,  early  experiments  with  various  test  masses  showed  that  the 
moving  mass  could  not  be  seen;  i.e.,  the  measured  force  was 
dominated  by  diaphragm  stiffness,  even  though  the  diaphragm  had 
been  made  as  compliant  as  practical.  To  resolve  this  condition,  a 
means  to  measure  pressure  in  front  of  the  piston,  inside  the  sample 
cavity,  had  to  be  developed.  Relocation  of  the  pressure  sensor 
permitted  a  stiffer  diaphragm  to  be  employed,  and  thereby  reduce  the 
effects  of  diaphragm  flexure  on  the  acoustic  pressure  in  front  of  the 
piston.  The  Initial  diaphragms  were  made  of  0.001  in.  and  0.002  in. 
thick  Mylar,  and  the  final  diaphragm  was  constructed  from  0.003  in. 
thick  beryllium  copper. 
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Figure  9.  Direct  Method  Equipment  Set-up. 


Several  alternative  methods  to  measure  the  pressure  inside  the 
sample  cavity  were  tested.  Fig.  9  shows  the  equipment  set-up.  For 
these  tests,  the  sample  cavity  was  filled  with  various  amounts  of 
Fluorinert  FC-72  non-polar  electronic  liquid  (3M  Corporation,  St.  Paul, 
Mn.).  Fluorinert  was  used  because  it  is  non-conductive  and  it  permits 
the  selection  of  various  viscosities  that  will  be  used  in  later  (post 
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thesis)  experimentation.  First,  a  variety  of  Polyvinylidiene  Fluoride 
(PVF2  or  PVDF)  films,  ranging  from  0.001  in.  to  0.01  in.  thick 
(Pennwalt  Corporation,  Folsom,  Ca.),  were  tried.  Each  sensor  was 
bonded  to  the  diaphragm.  In  every  case,  the  resulting  voltage 
response  from  the  pressure  sensor  was  contaminated  by  resonances  of 
flexural  waves  in  the  diaphragm.  Next,  a  similar  series  of  trials  was 
conducted  using  PVDF  sensors  suspended  in  the  fluid  column  in  the 
vicinity  of  the  piston  drive.  Again,  flexural  waves  dominated  the 
response,  in  this  case  propagating  in  the  sensors  themselves.  A  0.001 
in.  thick  electreted  PTFE  Teflon  sensor  was  tested,  with  no 
improvement  over  the  PVDF  series.  Finally,  a  small  cylindrical  tube 
piezoceramic  hydrophone  6  mm  in  diameter  was  tested.  The  small 
hydrophone  was  selected  because  its  lowest  frequency  resonance  is 
well  above  the  highest  frequency  of  interest  in  this  experiment.  While 
the  hydrophone  lacked  the  sensitivity  of  the  PVDF  sensors,  the 
measured  pressure  spectrum  over  the  range  100  to  2000  Hz  was 
substantially  improved. 

Figure  10,  raw  data  output  from  the  HP  4194  for  80  ml  added 
mass  of  Fluorinert,  is  an  example  of  the  voltage  response  of  the 
accelerometer  (lower  trace)  and  the  hydrophone  (upper  trace)  as  a 
function  of  frequency,  plotted  on  logarithmic  scales.  The  verticle 
scale  for  each  trace  spans  five  decades;  the  minimum  and  maximum 
values  are  listed.  The  verticle  scale  markings  are  for  the 
accelerometer  (A)  trace.  The  frequency  scale  ranges  from  30  to  3000 
Hz.  The  signal  to  noise  ratio  of  the  accelerometer  output  (lower  trace) 
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Figure  10.  Representative  Response  of  Accelerometer  ( "A:  R”  .  lower  trace  ) 
and  Hydrophone  ( "B:  T’ ,  upper  trace  )  . 
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is  very  high,  while  the  hydrophone  output  (upper  trace)  signal  to  noise 
ratio  is  poor  at  the  lower  frequencies.  The  accelerometer  output  is 
roughly  proportional  to  f2  for  frequencies  above  200  Hz,  indicating  a 
constant-displacement  drive  in  this  region.  Below  200  Hz  the 
acceleration  is  largely  constant.  For  the  most  part,  the  hydrophone 
output  signal  is  nearly  proportional  to  the  accelerometer  output,  the 
signature  of  a  mass- like  load  presented  to  the  driver.  The  peak  in  the 
hydrophone  response  at  about  2000  Hz  is  a  broad,  contaminating 
resonance  of  unknown  origin.  As  the  use  of  a  hydrophone  is  not 
compatible  with  the  closed  apparatus,  it  is  regarded  as  an  interim 
solution. 

Because  of  the  low  signal  to  noise  ratio  observed  at  low 
frequencies  (10-100  Hz),  considerable  effort  went  into  eliminating 
noise  sources.  Though  a  high  driving  voltage  (^  100  V)  was  desired, 
the  available  high  voltage  amplifier  was  found  to  introduce  excessive 
noise  into  its  output.  A  lower  voltage  output  amplifier,  Techron  7520 
Power  Supply  Amplifier  (Crown  International,  Inc.,  Elkhart,  In.),  was 
selected  for  its  low  noise  output.  The  Techron  amplifier  is  limited  to 
approximately  25  V  rms  output  voltage.  Pre-amplifiers  were  operated 
on  battery  power  to  avoid  60  Hz  harmonic  corruption  in  the  pressure 
and  acceleration  signals.  An  exhaustive  effort  to  eliminate  ground 
loops  among  various  equipment  (in  Fig.  9)  proved  to  be  most 
beneficial. 
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3.  Drive  Linearity 


Prior  to  collecting  calibration  data  on  the  direct  method 
apparatus,  a  characterization  of  the  drive  linearity  was  performed. 
Table  3.1  lists  the  accelerometer  and  force  gauge  outputs  of  the 
impedance  head  for  applied  voltages  to  the  piezoceramic  stack  of  2V 
and  20V.  The  accelerometer  and  force  gauge  voltage  responses  were 
collected  by  the  gain-phase  analyzer.  Data  at  each  frequency  were  also 
evaluated  for  the  presence  of  non-linearities  using  an  HP  3580A 
Spectrum  Analyzer.  It  was  observed  that  the  amplitudes  of  all 


TABLE  3.1.  DRIVE  LINEARITY 

Freq.  (  Hz.) 

Acceleration  (  pV) 

2  V 

20  V 

200 

16.5 

162 

498 

65.5 

651 

1005 

246 

2457 

2005 

1023 

10230 

Freq.  (  Hz.) 

Force  (mV) 

2  V 

20  V 

5.94 

59.4 

498 

5.89 

58.9 

1005 

5.67 

56.7 

2005 

4.74 

47.4 
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harmonics  were  at  least  30  dB  below  that  of  the  fundamental  driving 
frequency.  As  a  result,  it  was  concluded  that  the  drive  is  linear. 


The  objective  of  this  final  experiment  was  to  obtain  a  set  of 
calibration  data  for  the  apparatus  which  can  be  subsequently  used  to 
determine  the  unknown  moving  mass  of  a  fluid  contained  in  a  rigid 
porous  solid.  An  evaluation  of  the  calibration  data  set  was  made 
possible  by  a  comparison  with  results  for  known  moving  masses  of 
fluid.  For  a  nearly  constant  drive  voltage  of  25  V,  calibration  pressure 
and  acceleration  measurements  were  recorded  as  a  function  of 
frequency  using  the  equipment  in  Fig.  9.  Fluorinert  was  added  in 
increments  of  10  ml  ,  from  50  to  90  ml  (the  specific  gravity  of 
Fluorinert  FC-72  is  1.68,  the  sound  speed  is  525  m/s).  As  with 
previous  experiments  using  water,  the  viscous  penetration  depth  of 
FC-72  is  much  smaller  than  the  chamber  sleeve  (pore)  radius;  e.g.,  at 
100  Hz,  5  fc-72  =  10°  H.  Plots  were  obtained  for  pressure  and 
acceleration  for  the  band  20  to  2000  Hz  at  each  10  ml  increment  for 
qualitative  analysis,  e.g..  Fig.  10.  Because  a  broad  resonance  around 
2000  Hz  dominated  the  upper  end  of  the  frequency  range,  concurrent 
recordings  of  pressure  and  acceleration  were  taken  at  regular 
frequency  intervals  only  over  the  range  100  to  1000  Hz.  The  data 
were  collected  using  long  time  integration  and  16  averages  (features  of 
the  HP  4194),  and  the  resulting  calibration  data  set  was  stored  in  table 
format  (a  convenient  function  of  the  HP  4194). 


Drive  voltage  uncertainty  was  less  than  2%  and  added  volume 
uncertainty  was  less  than  4%.  Most  of  the  added  mass  error  was 
attributed  to  evaporation  of  Fluortnert  FC-72  at  room  temperature. 
Since  extracted  mass  is  a  quotient  of  the  two  experimentally  measured 
values,  per  Eq.  (2.17),  each  known  to  better  than  4%,  overall 
experimental  error  is  less  than  6%. 

Using  60  and  80  ml  data  as  the  calibration  control  set, 
extraction  of  the  moving  mass  for  50,  70,  and  90  ml  were 
accomplished  at  the  selected  frequencies  from  the  following 
interpolation/extrapolation  formula: 


Volume  extracted  ( co  )  =  60  +  (80  -  60) 

where,  for  a  given  mass  loading  and  frequency:  Vp  iS  the  output  voltage 
of  the  pressure  sensing  hydrophone  and  Va  is  the  output  voltage  of  the 
accelerometer.  The  results  of  the  extraction  of  the  volume  of  moving 
mass.  Fig.  1 1 ,  are  well  within  experimental  error  over  most  of  the 
decade  in  frequency,  100  to  1000  Hz,  and  constitute  a  successful 
calibration  of  the  direct  method  apparatus.  Significantly,  the  ability  to 
discern  the  volume  of  unknown  added  mass,  or  the  resolution  in  ml,  is 
5%  or  about  4  ml. 
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VOLUME  of 


RESULTS  :  EXTRACTED  MASS  vs.  FREQUENCY 


FREQUENCY  ( Hz ) 


Figure  11.  Calibration  Results  for  the  Direct  Method. 


IV.  INTERPRETATION  OF  RESULTS 


A.  COMPARISON  OF  TECHNIQUES 

It  is  important  to  compare  the  preliminary  calibration  results  from 
the  two  methods  and  identify  the  more  promising  technique  for  future 
investigations  with  porous  media  of  unknown  complex  mass.  The 
impedance  method  is  procedurally  very  easy  to  perform  with  the  HP 
4194  Impedance/Gain-Phase  Analyzer.  However,  because  the 
extraction  of  moving  mass  is  indirect,  the  impedance  method  is  a 
technique  which  requires  much  finesse;  a  precise  knowledge  of  the 
properties  of  the  transducer  is  essential.  What  the  direct  method 
lacks  in  experimental  ease  is  made  up  in  its  direct  extraction  of 
complex  mass. 

Apparatus  prototypes  for  both  techniques  were  developed  and 
employed  successfully.  The  preliminary  results  of  the  impedance 
method,  represented  by  Fig.  4,  are  comparable  to  those  of  the  direct 
method,  Fig.  11.  Five  per  cent  uncertainty  in  extracted  mass  over 
most  of  the  frequency  band  100  to  1000  Hz  was  established  in  both 
techniques.  However,  because  the  impedance  method  is  plagued  by  a 
systematic  deviation  in  extracted  mass  at  higher  frequencies,  which 
increases  in  magnitude  with  increasing  added  mass,  it  is  doubtful  that 
moving  masses  greater  than  50  g  can  be  extracted  using  this 
technique  with  the  present  moving  coil  driver.  The  direct  method 
exhibits  no  such  error.  Uncertainty  in  applying  the  direct  method  can 
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further  be  reduced  by  refinement  of  the  design  of  the  piston  and 
surrounding  diaphragm  to  reduce  the  effect  of  flexural  waves. 


B.  THE  PREFERRED  METHOD 

Based  on  the  preliminary  results  of  both  experimental  techniques 
to  date,  the  use  of  the  direct  method  for  further  Investigation  is 
recommended.  While  there  is  room  for  improvement  in  the  apparatus 
of  both  techniques,  the  impedance  method  is  limited  by  a  driver  that 
appears  unsuitable  for  extraction  of  moving  masses  greater  than  50  g. 
Conversely,  it  has  been  demonstrated  that  it  is  possible  to  extract  a 
moving  mass  of  approximately  the  value  expected  for  fluid-filled 
porous  solid  samples  (~  100  g)  with  a  resolution  of  several  percent 
using  the  direct  method. 
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V.  CONCLUSIONS 


An  experimental  apparatus  was  developed  by  which  the  complex 
mass  of  a  viscous  fluid  contained  in  a  rigid  porous  solid  can  be  directly 
extracted  from  acoustic  pressure  and  acceleration  measurements. 

Several  preliminary  experiments  were  conducted  to  verify  the 
ability  to  extract  the  value  of  a  known  mass  of  fluid.  The  results  for 
the  extracted  mass  were  accurate  to  within  five  per  cent  of  known 
values  over  most  of  the  frequency  range  100  to  1000  Hz. 

A  comparison  of  the  preliminary  results  of  the  two  techniques  of 
obtaining  complex  mass  demonstrated  that  the  direct  method  was 
most  suitable  for  use  in  further  work  investigating  porous  solid 
geometries. 
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APPENDIX  A  OPEN  SYSTEM  THEORY 


Calibration  of  both  the  direct  and  impedance  methods  was 
accomplished  using  the  open  system,  single  transducer/test  chamber 
set-up  of  Fig.  9  (direct  method)  and  Fig.  12  (impedance  method). 
The  effect  of  the  radiation  load  presented  to  the  open  end  on  the  load 
presented  to  the  driver  needs  to  be  examined.  Ignoring  the  radiation 
load  for  the  moment,  in  the  long  wavelength  limit,  the  impedance 
presented  to  the  face  of  the  driver  by  the  liquid  in  the  chamber  alone 
is  (Grant,  1988,  pp.  8-10), 

^radomvEji  ~  ^rad  liq  -j®Mr  j  |^_  JJ 

where  Mr  is  the  mass  of  the  liquid.  Since  the  particle  velocity  is  the 
same  at  the  top  and  bottom  of  the  liquid  in  the  low  frequency  limit, 
the  impedance  seen  by  the  driver  is  Just  the  sum  of  the  above,  Eq. 
(A.  1),  and  that  of  the  radiated  sound  wave  in  air: 

^rad drive*  —  ZraduQ  ^radM*  #  (A.2) 

Substituting  Eq.  (A.  1)  for  Zraduo  and  the  appropriate  expression  for 
Zrad*i*  ,  assuming  the  worst  case  of  ka  =  1,  into  the  right  side  of  Eq. 
(A.2)  gives: 

Zrad«,rvE*  =  jwMr  +  PcS  ,  (A.  3) 
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where  S  is  the  area  of  the  sample  at  the  fluid/air  boundary.  Since 


M.  «i 

l)“MJ  (A.4) 

for  the  stated  experimental  conditions  and  dimensions,  Eq.  (A.3)  may 
be  approximated  as: 

=  COMf  =  [ZuuIljqI  .  (A.  5) 

Thus,  the  theory  in  Chapter  II  does  not  need  to  be  modified  for  the 
open  system. 


Figure  12.  Single  Transducer /Test  Chamber  Open  System. 
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APPENDIX  B  TRANDUCER  CALIBRATION  PROGRAM 


CALDUCER. 

*«*****»»»»*»****#*★*#******#♦*►*»*»****##*♦***#***#*###♦*♦*►*##.*  *.* 

DIRECT  MEASUREMENT  METHOD  -  COMPLEX.  MASS  DENSITY 
CALIBRATION  PROGRAM 

ir**»ifr'**»»**ir*^*,iMp#*i**fc#**ir***##ifri*iMMM**#i»,'*iMMfr*#*#*«'*#****4f*»*##*',fr 


CONTROL  !  ,12.0 
PRINTER  IS  CRT 
PRINT  USING  "?’" 
PRINT  USING  " / ' " 


!  So f treys  ON 


w:tten  by  R.  A.  MIRICK 


Lest  update  DO  NOV  SS 


DISR  "Ensure  rawer  ON  to  Analyzer,  DAC  ,  and  ail  ancillary  equipr 
WAIT  5 . 

PRINT 

DI5F  "Press  CONTINUE  when  ready  ..." 

PAUSE 


►  ****■»•*•»********♦***. 


1  INITIAL  CONDITIONS 

nnrrnu  0  *  r-  r-  , 

yr  i  i u'v  i 

INTEGER  I  ,Np t s 

D I M  Prrq_h z <  ZOO  / .Acce  1 7ml.  s  *'  ZOO  ) , Acce  1  ~aeg >'  ZOO  ' ,  Acce 1 6m>  s  (  ZOO  '• 
DIM  Acce  15 deg ( 200 /  .ForceSmk s »  ZOO  ? , Force Sdeg (  ZOO  )  .Force  1  Omks  <  ZOO 
DIM  For ce "i Odeg >  ZOO  > , Re f_mdb v(  ZOO  ) , Bn i nus_a~negdb '  ZOO  < 

DIM  Dniinu5_sSnagdb(  ZOO  /  ,  Accei_7nagv<  ZOO  >  ,Accei__8magv(  ZOO’  ' 

DIM  Bmx nus_a5nagdb '  ZOO  >  ,Bminus_e  1  Omgdbv  ZOO  >  ,Force_5n.agv(  ZOO  ) 

DIM  For ce_ i  Ong v (  ZOO ’’ 

MASS  STORAGE  IS  "  ,700,'“ 

DEG 

0 s c_ 1 e ve i “ . 5  i  Volts  rns 

ASSIGN  GAneiyzer  TO  7i7  1  Analyzer  “  HP415ZA 


ASSIGN  ®Dsc  TO  70S 

CLEAR  BAnaiyzer 

CLEAR  (Suae 


)AC  “  HF341 


STATUS  OAnaiyzer; Srq_s  t a toy t e 
IF  Src_statbyte“S4  OR  5rq_statbyte>64  THEN 
uisr  iNrOi  trtrur  .  ,  • 

GOTO  i S i 0 
END  IF 


trror  eust: 


TERM I NATE 
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a  i  <o 

520 

r—  —  /I 

201 a 
540 
550 
55  0 
570 
580 
5S0 
600 
6  ( 0 
520 

5  30 
640 
650 

6  50 
670 
680 

6  30 

7  00 
7  i  0 
720 
730 
740 
750 
750 
770 
780 
730 

m  /-»  /-» 

aw 
6  1  0 
620 
530 
840 
550 
860 
670 
58G 
6  30 
300 
3  c  0 
S20 
530 
340 
350 
360 
570 
350 
330 

f  000 


j*ir»***»#'***»***‘i**»**ir***>*»tr***'**‘»«**##*#,»*»'if#**#*»**##**#***»'»*** 

!  CALIBRATION  PROCEDURE 

lirir*****#»irir*iri»ir*****ir**iMr**ir*****^****************^if***###«****** 

OUTPUT  ©Analyser:  "CL''  ;Gsc_level ;  “EN" 

OUTPUT  ©Anaivzer :  “N2T2V1  “  1  External  trigger,  measure  in  d8V , 

!  AVERAGE  ON 

OUTPUT  ©Analyser: "EX " 

WAIT  4 . 

Np t  s_8 

FOR  I“1£  TO  2008  STEP  ! 0  •  I*  Frequency  in  Hz 

Npts“Npts+ t 
Freq_K  hz“I  X i 000 . 

©Analyzer ; "FR" ;Freq_khz; “EN"  I  Spot  frequency  set  to  I  H 
„  !  Allow  settling  time  for  impedance  follower 

I  operational  amplifier 

BEEF  30 i i  .  , . 5 

'FREQUENCY  IS  "  ,1,"  Hz" 


OUTPU 
!  WAIT 


©Ana i y zer : " A6N2V t 
©Anaivzer : " EX " 


Function:  measure  REF  input 

Trigger  command 

Take  advantage  of  AVERAGE-ing 


©Ana iyzer ; Re f_indb v( Np t s ) 

Reference  input  is  "  ,Ref_mdbv<  Npts  /  , "  dBV" 


uicr 
WA IT  2 . 

OUTPU 
OUTPU 
WAIT 
ENTER 
01 SF 
WAIT  4. 

iQISP  "Observe  input  on  Qscope 
OUTPUT  ©Analyser  ••  "nENZv'i  " 

OUTPUT  @uac;"OrN" 

OUTPUT  ©Dac;"CLS7" 

WA IT  .01 

OUTPUT  ©Ana  1 v zer ; " EX " 

WA  I T 

©Ana  1 y zer : om inus_a7magdb ( Np  t s  )  ,  Acce 1 7 deg'  Np  t  s  > 
suae : "0FN7" 

©Dac : "CL Sc " 

G  i 

©Ana  1 v:er : "EX " 


Ensure  signal  is  not  clipped 
Function:  measure  TEST  input 
1  Open  ail  relays 

1  Ensure  reiav  closure 


ENTER 
OUTPU" 

OUTPUT 
WA  I T 
OUTPUT 
WA  I T 

ENTER  ©Ana i y zer ; 6m inus_aomagdb (Npts  )  ,  Acce 1 8 deg >  Np  t  s > 

OUTPUT  ©Dac;"GPN3" 
i  -  -  -  dBv  to  volts  i 

Acce i_7magv( Npts  >-1 G . " ( ( Bminus_a7magdb (Npts  )+Ref_indbv( Np t s  /  it  1 
Acce l_Smagv ( Np t s  /“ i 8 . “ ( ( Bmi nus_aSmagdb ( Np ts  )+Ref_indbv( Np t s  )  /XI 
1  -  -  -  Volts  to  mX  sec "2  ! 

Acce  17mf.s  >  Np  t  s  '/-Acce  l_7mag v(  Npts  >*303.5 
Acce  iSmf.s  (  Np  t  s  >"Acce  i_Bmagv(  Np  t  s  /*324.7 
i-  -  -  NOTE  :  THIS  CONVERSION  APPLIES  ONLY  TO  IH#  3287  a  3315  -  -  - 

Ca l_di f  f “Acce 1 7mks ( Npts  >XAcce ISmks < Np t  s > 

DI3F  "  Acceleration  ratio  is  ",Cai_diff 


40 


6) 


4 


% 


(  'a  ( 

1  020 
r  030 
[040 
1  050 
1  05 0 
1  070 
1  080 
1  050 
It  00 
i  1  1  0 
1  !  20 
1  i  30 
1  1  40 
t  1  50 
m  S3 

1  f  70 
It  50 
1  f  50 
1  230 
1  2 ' 0 
•220 
1  23 G 
f  240 
(  250 
*.250 
i  270 
I  280 
1  250 
l  300 
(310 
1323 
l  333 
i  340 
1  350 
1  353 
1  373 


wn I 7  2 . 

IF  C4l_di.fr >1.1  OR  Gai_diff-C.S  THEN  I  Accuracy  of  procedure 

0 ISP  "  WARNING...  Acceleration  ratio  13  "  ,Cal_diff 
WAIT  4 . 

INPUT  "Go  you  wish  to  try  again7  Y/N  ” ,Answ$ 

IF  AnsuS*"Y"  OR  Answ$“"yM  THEN  GOTO  780 
END  IF 

OUTPUT  @Dac ; “ 0PN3 M 
OUTPUT  ?Qac ? "CLSo” 

WAIT  .01 

OUTPUT  ryAna  I  veer  ;  "  EX” 

WAIT  *  . 


ENTER  £Ana  Ly  zer  :  Bn  inus_*Snagdb  (  Np  t  s  >  ,Force3deg(  Np  t  3  / 
OUTPUT  tyOac  ;  " urNS  ’* 

OUTPUT  @uac ; " CLS 1 3 " 

WAIT  .01 

OUTPUT  ^Analyzer EX " 

WAIT  I  . 


cim  (  err  vnnaiyzer  ;  on  mus_a  t  tamgao  v  ivp r. 3  /  7rorce  1  <aaegv  wp  t  3  / 

OUTPUT  80ac  5  H0PN 1  0*' 

-  -  -  dSV  to  Volts  -  -  -  1 

Force_Snagv «  Np  t  s  > » 1  0 .  (  ■'  Bn inus__a9nagdb  (Nets  >+Re f_mdbv(  Np  t  s  )  / /  1 
r  orce_  1  Ongv*.  Np  t  3  /*  1  3.  A  «  (  Sn  1  nua_a  1  0ngdb  (  Np  t  3  /+Re  f_i  ndbvi  f’io  t  3  )  ) /  1 

-  -  -  volts  to  Newtons  -  -  » 

For ce Snr s ( Np  1 3  >-Force_Snagv<  Np  1 3  >*2 . 525 
Force  1 3nr 3  v Np  1 3  /■rorce^i 0ngv<  Np  t  s >*Z . 538 
-  -  NOTE  1  THIS  CONVERSION  APPLIES  ONLY  TO 
Fre q_hz<  Np  t  5  /*  I 


TO  IH4  525 


a  <  <1  0010 


(VC.  A  I  1 

STATUS  8Ana i y zer ; Srq_s t a t b v t e 
IF  Srq_s tatby t e-54  OR  Srq_s tatay t e>54  THEN 
ulSr  “INPUT  ERROR..." 

GOTO  1 5 1  0 
END  IF 


1  Error  exists' 
I  TERMINATE 


* 


4 1 


L 


cs>  cs> 


fc 


0 


483 

}  ir  •# 

*  *  *  *  * 

413 

t 

OUT 

PUT 

428 

J  *  •* 

*  ♦#  -4 

♦  *  ir  ir  * 

1430 

INF 

u7 

Enter 

t  he 

443 

CRE 

nTE 

GOAT 

NameS 

453 

nSS 

IGN 

©Fiie 

_name 

453 

Fur 

I  »  < 

7Q  Npta 

473 

OU 

TPUT 

©F  i  ie 

480 

QU 

7PU7 

©Fiie 

430 

OU 

7rU7 

@F  lie 

533 

QU 

7PU7 

©F  i  ie 

5 1  3 

QU 

7PIJ7 

©Fiie 

5Z3 

uu 

TPUT 

©F  i  ie 

533 

OU 

TPUT 

©Fiie 

34  8 

OU 

7PU7 

®F  i  Ie 

1  553 

OUTPUT 

©Fiie 

553 

NEX 

i  L 

573 

nS  5 

IGN 

©File 

_nane 

►  #**** 
.NameS 


■iririri*ifrir*'*ifr*#*ifrir*#****'*irir*'* 


me  ;rreq_nz<  i  t 
©F i i«_naroe ; ncce I7nk s < I 
®F 1  ie_name;nccei7deg(  I 
©Fiie  name  s ncce  13mka  I 


i  ooio 

«  r»  /n i 

!  33U 

1  503 
f  5  t  3 
1  5Z3 
i  533 
I  5«2 


.me  ; ncceiaaegv  i  » 
@F i Le_name ; Forceumk a <  I > 
me ?Farce3aeg(  I > 
.me  ;  Force  I  3mk  a  (  I  > 
:me  ;  Force  i  0deq(  I  > 


t  u  * 

»*ir*ir*i**ir*-*ir‘iririr'*-irieir-eiri**ifr****i*ir*ir'**ifr*iririr'«r,< 

TERM INATE 

PRINT  USING  “©“ 

PR  INT  US  INS  "  f  t  !  " 

PRINT  "Program  conclusion...'1 
F  R I N  T 
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